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Laser-Pumped  Coherent  X-Ray  FEL 


P.  Sprangle,  B.  Hafizi+,  and  J.R.  Penano 


Plasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  DC  20375 
+Icarus  Research,  Inc.,  PO  Box  30780,  Bethesda,  MD  20824-0780 


Abstract 

In  a  laser-pumped  x-ray  free  electron  laser  (FEL)  an  intense  laser  field  replaces  the 
magnetic  wiggler  field  of  a  conventional  FEL.  Depending  on  the  intensity  and  quality  of 
both  the  electron  beam  and  pump  laser,  the  Thomson  backscattered  radiation  can  be 
coherently  amplified.  In  a  conventional  FEL  the  generation  of  x-rays  requires  electron 
beam  energies  in  the  multi-GeV  range.  In  a  laser-pumped  x-ray  FEL  electron  beam 
energies  in  the  multi-MeV  range  would  be  sufficient.  To  generate  coherent  x-rays  with 
this  mechanism  a  number  of  physics  and  technology  issues  must  be  addressed.  Foremost 
among  these  are  the  stringent  requirements  placed  on  the  electron  beam  quality  and 
brightness  as  well  as  the  pump  laser  intensity  and  pulse  energy.  The  seed  radiation  for 
the  laser-pumped  FEL  is  the  laser-induced  spontaneous  radiation.  The  evolution  of 
incoherent  radiation  into  coherent  radiation  as  well  as  the  power  gain  lengths  associated 
with  the  coherent  x-rays  are  analyzed  and  discussed.  There  is  excellent  agreement 
between  our  analytical  results  and  GENESIS  simulations  for  the  radiated  power,  gain 
length,  conversion  efficiency,  line-width  and  saturation  length.  These  issues,  as  well  as 
others,  necessary  to  achieve  coherent  amplified  x-rays  in  a  laser-pumped  FEL  are 
discussed. 
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I.  Introduction 


The  free-electron  laser  (FEL)  can,  in  principle,  generate  coherent,  polarized,  short 
pulses  of  x-rays  for  numerous  applications  in  research.  There  are  a  number  of  large-scale 
electron  accelerator  facilities  throughout  the  world  which  will  be  used  tor  x-ray 
generation  using  a  conventional  FEL  configuration  [1-4].  In  a  conventional  FEL  the 
electron  beam  propagates  through  a  static,  periodic  magnetic  field  (wiggler)  which  results 
in  stimulated  emission  [5-20].  Generation  of  x-rays  at  these  facilities  typically  requires 
electron  beam  energies  in  the  multi-GeV  range  with  peak  currents  in  the  multi-kA  range, 
and  wiggler  lengths  of  many  tens  of  meters.  An  x-ray  FEL  amplifier  c£n  be  operated  in 
the  self-amplified  regime,  eliminating  the  need  for  a  coherent  input  x-rtiy  source  [21-26]. 
In  this  case  the  FEL  seed  radiation  is  provided  by  spontaneous  incoherent  emission  in  the 
wiggler. 

The  wiggler  field  in  the  FEL  can  be  replaced  with  an  electromagnetic  wave  such 
as  an  intense  laser  field.  Early  analysis  of  stimulated  emission  from  rel  itivistic  electrons 
interacting  with  an  electromagnetic  pump  was  presented  and  discussed  in  [27],  This 
analysis  was  limited  to  the  low-gain,  thermal  beam  regime.  In  this  regime  the  power  gain 
lengths  are  extremely  long  making  the  concept  impractical.  The  high-gain  regime  of  the 
electromagnetically  pumped  FEL  was  first  analyzed  and  discussed  in  [28].  In  this  regime 
the  power  gain  lengths  can  be  very  short.  However,  the  requirements  on  the  electron 
beam  quality  and  the  pump  laser  power  are  demanding,  particularly  for  x-ray  generation. 
Since  these  early  studies  there  have  been  a  number  of  papers  that  have  considered 
employing  electromagnetic  pumps  in  FELs  [29-32], 
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In  this  paper  we  analyze  and  discuss  a  laser  pumped  FEL  amplifier  operating  in 
the  x-ray  regime.  The  analysis  considers  i)  electron  beam  thermal  effects,  ii)  off-axis 
propagation  and  iii)  the  transition  from  incoherent  (spontaneous)  to  coherent  x-rays.  The 
power  gain  length  and  the  conversion  efficiency  are  determined  as  functions  of  the 
electron  beam  energy  spread.  The  radiation  power  as  a  function  of  interaction  distance  is 
obtained  in  both  the  incoherent  and  coherent  regimes.  The  coherent  power  is  emitted  into 
a  solid  angle  which  is  typically  much  greater  than  the  solid  angle  associated  with 
diffraction.  For  electron  beams  of  sufficiently  high  quality,  with  energies  of  7MeV  and 
peak  currents  of  500  A  ,  we  find  that  coherent  x-rays  at  15  A  can  be  generated  with 
power  gain  lengths  of  500// m,  saturation  lengths  of  0.4cm  and  conversion  efficiencies 
of  ~  0.01  %.  To  achieve  these  values  the  fractional  electron  beam  energy  spread  must  be 
<  0.01  %  .  The  pump  laser  for  this  example  has  a  wavelength  of  1pm,  pulse  duration  of 
27psec  and  power  of  50TW  .  To  compare  our  results  with  simulations  we  use  the 
GENESIS  FEL  code  [33]  and  find  good  agreement  with  our  analytical  results.  We  also 
compare  our  results  and  GENESIS  simulations  for  the  LCLS  wiggler-based  FEL 
operating  at  15 A. 

II.  High-Gain  Regime 

The  laser  pumped  FEL  is  shown  schematically  in  Fig.  1.  The  pump  laser  is  taken 
to  be  circularly  polarized,  with  normalized  vector  potential 

a„(r,/)  =  an  (cos  (k^z  +  C0ot)ex  +  sin  (k„z  +  aot)ey) 

=  (a0/y[2)exp(-i(k0z  +  coj))eL  +  c.c. , 

where  ex  =  (e^  +/ev)/V 2  is  a  unit  transverse  vector,  Atl  =  2kcI(00  is  the  pump 
wavelength,  k(1  =  (0o/c  is  the  wavenumber  and  a0  =  q  A0 /me2  is  the  normalized 
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amplitude.  The  pump  laser  power  is  P0  =  (m2 c5 / q2)(n<J0 /  A2) a],  where 

m1  cs  / q2  =  8.75  GW ,  oo  =  n  r]  /  2  is  the  cross  sectional  area  for  a  Gaussian  transverse 

profile  and  ra  is  the  laser  spot  size. 


The  x  ray  radiation  is  given  by  the  normalized  vector  potential 
a(r,0  =  [a(r,t)/^)e\p(i(kz  z  -  cot))e±  +  c.c.. 


(2) 


where  A  =  2nd co  is  the  x-ray  wavelength  and  k.  is  the  complex  axial  wavenumber. 
a)  3-D  Thermal  Beam  Dispersion  Relation 

Thermal  effects  associated  with  the  electron  beam  play  a  critical  role  in  the  FEL 
interaction.  The  FEL  dispersion  relation  including  thermal  effects  is  [13,15,16], 


drFM) 


(3) 


[kz+k<l-(G)-oi>)/\z  +juk[f  ' 


where  k±  is  the  transverse  wavenumber,  /  is  the  filling  factor,  i.e.,  ratio  of  electron 
beam  to  radiation  beam  areas,  v  -  o^r^/Ac2  =  Nhre/lh  =  7fr[A]/l 7,000  is  Budker’s 
parameter,  a \  =  (4 nq2nhl m)'12  is  the  electron  beam  plasma  frequency,  Ih  is  the  beam 
current,  rt  =  q'  I  me2  is  the  classical  electron  radius,  ib  is  the  electron  bunch  length,  Nh  is 
the  number  of  electrons  in  a  bunch,  rh  is  the  electron  beam  radius  and  Fn(y)  is  the 
electron  distribution  function.  The  filling  factor  is  a  function  of  the  interaction  distance. 
In  Eq.(3)  ju  =  y[0a20  /(2 y^oalc)  is  a  correction  term  that  arises  from  the  transverse 
electron  motion  in  the  field  of  the  pump  laser.  The  resonant  frequency  is  a  function  of 
kL ,  and  for  a20  «1  is  given  by 


(^1 )  —  0  k±)  1 


(4) 
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where  the  resonant  frequency  for  on-axis  ( kL  =  0)  propagation  is 

(Ot jo  =  4  yjpCO,,  =  4  yf}  con  /(I  +  al ) .  The  coherent  radiation  is  emitted  along  the  z-axis 

inside  a  narrow  cone  with  opening  angle  0k  =  k± / k . .  The  range  of  allowed  k±‘ s  ,  i.e., 

emission  solid  angle,  is  important  in  determining  the  incoherent  and  coherent  x-ray 
power  and  is  discussed  in  Sec  He. 
i)  Cold  Beam 

For  a  cold  beam  the  electron  distribution  function  is  F(t{y)  =  8{y-yn)  and  the 


dispersion  relation  is  given  by  Dfel(k,co)  =  0  where 

+(r„/V3)’.  (5) 

For  a  cold  beam,  with  kL  =  0  and  0)  =  0)m ,  we  find  k.  =  co/c  +  Ak ,  where 
Ak  =  (1/V3  -  /) T  / 2 .  The  power  growth  rate  as  a  function  of  (0  and  kL  is 

Tt(0),kj  =  rs„(l  -  {(O-  ajR(k±))2/Aco2),  (6) 

where  the  peak  growth  rate  is  1“  =  (5.07/ y,)(f  Va20  A.J  j'3 ,  the  power  gain  length 

is  L  =  1  /  r  ,  and  the  line  width  associated  with  the  power  growth  rate  is 

A co/ O)R0  =  (A()/  Ls(l)/  2k  .  As  an  example,  the  power  gain  length  for  x-rays  at  X  =  15  A 

is  shown  in  Fig.  2  for  a  cold  electron  beam  as  a  function  of  beam  current.  The  parameters 
for  this  plot  are  listed  in  Table  I. 
ii)  Thermal  Beam 

For  a  thermal  electron  beam  with  distribution  function 
F„{y)  =  (4ndyy  1  exp(-  (y  -  yo)2  /Sy2)  the  dispersion  relation  for  kL  =  0  is 


Dfel(k,o))= 


kz~ 


O) 


c2k 


lay, 


ro  y 


K-~ 

C 


1  + 


(Q)-Q)r o) 


fickl N 


2  T&to, 


R  0 


(O, 


RO  J 


5 


(7) 


Ak  =  — 
An 


\3 


dx  n  1/2  jcexp(-jc2) 


V  Yo_  r _ _ 

[Sj  Srlbk  -2ko{w-(oR0)/(oR0  +Ak0{5ylr0)x' 


where  dyt  jo  is  the  fractional  energy  spread.  The  dispersion  relation  can  be  written  in 
the  form 


#  =  -«,[i  +  !f+f„}z(;  +  #„)],  (8) 

where  £  =  -{yj dy)Ak! 4k0,  £0  =  {yJdy){co- (oR0)/2(oR0,  Ak  =  Akr  -zTg/2, 


P„  =2.4x10  5  (\  r^)3  {yo  /Syf ,  Z(£)  =  it~xn  J<7xexp(-.r  )/(*-£)  is  the  plasma 

dispersion  function  and  gain  occurs  when  the  imaginary  part  of  t;  is  positive.  In  the  cold 
beam  limit  +  £,|»1  the  dispersion  relation  reduces  to  the  usual  cubic  equation  with 
power  growth  rate  given  by  Eq.(6)  for  k±  =  0  and  =  0. 

The  thermal  dispersion  relation  can  be  analyzed  in  various  limits. 

In  the  thermal  beam  limit  <  1  the  dispersion  relation,  reduces  to  [15] 

£  =  -Po( 1  +  />/*(£  +  &)«p(-(£  +  O2)),  (9) 

where  Z(|£|  <  1)  =  zV^rexp(-£2) .  For  1>|£,|»|£|  the  imaginary  part  of  £  is 
=  -nU2po  £()  exp(— £2 )  where  po  «1 .  The  maximum  growth  rate  occurs  at 
£ 0  =  —  1  / -s/2  and  is  given  by  [15,27] 

rg/rgo  =  9.1x10 -*(A0rgo)2(rt,/#r)2,  (io) 

where  Ts.,  =  (5.0 1  /  yo)(f  v /(rf  A0)  J'  is  the  cold  beam  power  growth  rate  and  the 

thermal  growth  rate  is  inversely  proportional  to  the  square  of  the  energy  spread.  The 
power  growth  rate  at  resonance  (£,  =  0)  in  the  thermal  beam  limit  is 
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5.2X10-' 

i  +  1.8x10-’  aj„?(roiSrt' 


do 


r  /r  = 

J?  K» 


In  the  extreme  thermal  limit  |£|  «  1  the  power  growth  rate  at  resonance  is  given  by 

r,/r„  =  5.2xio-'a„r,„)>(r„/<%')5.  (i2) 

O 

Figure  3  plots  the  normalized  growth  rate  for  x-rays  at  X  =  15  A  as  a  function  of 
relative  electron  beam  energy  spread  and  detuning  for  the  parameters  listed  in  Table  I. 
Figure  3  shows  that  as  the  energy  spread  of  the  beam  is  increased  up  to  Syl  yo  -5x10  4 , 

the  FEL  interaction  can  be  detuned  to  increase  the  growth  rate  relative  to  the  resonant 
growth  rate.  For  a  given  energy  spread,  the  optimal  detuning,  i.e.,  maximum  growth  rate, 
occurs  when  the  difference  between  the  beam  velocity  and  the  phase  velocity  of  the  wave 
is  equal  to  the  thermal  velocity  spread  of  the  beam.  For  dyl  yo  >  5x  10”4  ,  the  growth  rate 
is  vanishingly  small  regardless  of  detuning. 
b)  X-  Ray  Conversion  Efficiency 

The  saturated  coherent  power  is  Pcohsa,  =  7]  Nh  y0m  c2  /  Th,  where  r]  =  Pcoh  sal  /  Ph  is 

the  conversion  efficiency  and  Ph  =  Nh  ynmc2  /rh  =v ynmc  I r  is  the  electron  beam 

power.  The  conversion  efficiency  in  the  cold  beam  limit  can  be  obtained  by  considering 
the  difference  between  the  electron  beam  energy  before  and  after  trapping  in  the 
ponderomotive  potential  [17].  The  efficiency  at  saturation  is 

//  =  (2 /7’0)(^7’,>/^Ao)(Ao  —  PPh)  >  where  (3.  is  the  normalized  axial  phase  velocity  of 
the  ponderomotive  wave.  From  the  dispersion  relation  the  phase  velocity  is  found  to  be 
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vph  lc  =  (co-coo)l c{kz  +kg )  =  J3,(>  -  R c(Ak)/k.  +  k\l(2k]).  The  conversion  efficiency 
at  saturation,  for  k±  =  0,  is  [17] 


7  =  0.023 (XJLj. 


(13) 


The  x-ray  conversion  efficiency  plotted  as  a  function  of  beam  current  is  shown  in  Fig.  4 
for  a  cold  electron  beam  along  with  results  from  GENESIS  simulations.  The  parameters 
for  this  plot  are  listed  in  Table  I.  There  is  good  agreement  between  theory  and 
simulations. 

c)  Validity  of  Classical  Description 

The  classical  description  is  valid  if  the  electron  momentum  recoil  is  somewhat 
less  than  the  electron  thermal  momentum  spread.  In  the  beam  frame  (irdicated  by  a 
prime  on  the  variables)  this  condition  is  tik'  «  m  Av'  and  in  the  laboratory  frame  it 
can  be  written  as  [34] 


(14) 


where  Ac  -  Intilmc  =0.02  A  is  the  Compton  wavelength, 

a y!y0  =  ^(l  +  ^)_1a  vz/c  is  the  fractional  energy  spread  and  Av.  is  the  axial  velocity 
spread. 

d)  Electron  Beam  Quality  Requirements 

For  high  gain  and  efficiency  the  electrons  must  remain  in  phase  with  the 
ponderomotive  (trapping)  wave.  As  a  result  the  interaction  is  sensitive  ::o  an  axial 
electron  velocity  spread.  A  spread  in  axial  electron  velocity  Av.  can  result  in  phase 
mixing  which  would  reduce  the  gain  and  efficiency.  The  electron  beam  can  be 
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considered  cold,  i.e.,  mono-energetic,  provided  Av_  Lgo  /  c  «A ,  which  can  be  written  in 


terms  of  the  fractional  energy  spread  Ay/ y«  A0  /  4  Lgo  =  10  7.  The  energy  spread  on 

the  beam  consists  of  several  contributions.  These  contributions  include:  i)  intrinsic 
energy  spread,  ii)  transverse  and  longitudinal  emittance,  iii)  space  charge,  iv)  pump  laser 
line  width,  and  v)  pump  laser  field  gradients.  The  overall  energy  spread  is 


Ay 

Y0 


A  y^ 

\Yoj 


+ 


A  y] 

\Y„j 


+ 


l.cmit 


Ay'' 

\Yoj 


+ 


'a£ 

\Yoj 


+ 


'Ay' 


space 

charge 


+ 


v  Yo 

y  iinewidlh 


'a£ 

Yo 


,  (15) 


pump 

grad 


where  (A y/yo)lmit  =  £ln  / 2r2  ~  2xl0"\  (e„  =  lmm-mrad,  rb  =  50 /tm), 

=  £n.J(TbEb)  ~  2.5x10"*  =25keV-psec  ,  rb  =  lOpsec,  and 


Eb  =  lOMeV ) ,  (A y/y  \p*x  =  |//v  (Ar/r.Vnp  =  SA/2A,  <  10"*.  Here,  £„isthe 

h  '  ■'charge  ' "  v  °  Miivwidlh  "  n 


normalized  transverse  emittance  and  e‘n ,  is  the  normalized  axial  emittance  [35].  The 

energy  spread  contribution  due  to  space  charge  leads  to  an  energy  shear  which  can  be 
eliminated  or  substantially  reduced  by  an  appropriate  tailoring  of  the  transverse  gradients 
of  the  pump  laser. 

The  electron  beam  brightness  Bn  =  2 Ih  l(7C2  e 2)  is  a  measure  of  beam  quality 
[35].  At  the  cathode  the  normalized  brightness  can  be  expressed  as 
Bn  =  Jcmc2  l(2/rkBTc)  where  is  the  current  density  at  the  cathode  and  Tc  is  the 

cathode  temperature.  For  a  photocathode  with  kBTc  =0.1eV  and  Jc  =  100  A/cm2  the 
brightness  is  Bn  ~  10s  A/(cm-rad)2 .  The  brightness  needed  in  a  laser  pumped  x-ray 

FEL  is  about  an  order  of  magnitude  higher,  i.e.,  ^lO1"  A/(cm-rad)2 .  An  axial  magnetic 
field  may  be  necessary  to  guide  the  electron  beam  though  the  interaction  region.  The 
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magnetic  field  required  fora  matched,  i.e.,  constant  radius,  electron  beam  is 
B[kG]  =  (4J/rb[cm])(v/yoy'2  ~  20kG. 
e)  Radiation  Solid  Angle 

The  transition  from  spontaneous  to  coherent  radiation  is  critically  dependent  on 
the  angular  distribution  of  the  waves.  Waves  with  finite  kj  have  a  propagation  angle 
6k  =  k±/kz  with  respect  to  the  z  axis.  The  peak  growth  rate  is  independent  of  k±  for 
waves  propagating  in  the  near-forward  direction  as  indicated  in  Eq.(6).  However,  as  k± 
increases  the  resonant  frequency  CO  -  (0R{kL)  decreases  as  indicated  schematically  in 
Fig. 5.  The  minimum  propagation  angle  is  0kmin  =  6 D  =  Xlnrs  where  0[}  is  the 
diffraction  angle,  rs  is  the  radial  dimension  of  the  radiation  beam  and  k  min  —  2 /  rs .  In 


general,  however,  k±  can  be  significantly  greater  than  k  mn .  From  the  power  gain 

expression  in  Eq.(6),  the  maximum  transverse  wavenumber,  for  gain  at  resonance 
( co=coR0 ),  is  given  by 


„  y/2 


k  -  0  k 

*±,max  uk, max 


& 


J 


K- 


(16) 


The  ratio  of  the  maximum  to  minimum  transverse  wavenumbers  is 


^i.max  / min  ~  3 (ZR  / L  )  where  ZR  =  7Tr2  /  X  is  the  Rayleigh  length  .  For  a  laser 

pumped  FEL,  ZR  »  Lgo.  In  the  case  of  an  optically  guided  FEL  amplifier  ZR  ~  Lgo . 

Other  processes,  such  as  the  electron  transverse  wiggle  and  betatron  oscillations,  can  also 
limit  the  range  of  transverse  wavenumbers. 
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The  solid  angle  associated  with  the  radiation  beam  is  ALlk  =  n where 

9kma  The  spontaneous  (incoherent)  radiation  is  directed  into  a  forward 

cone  with  angle  0incoh  ~  1  lyw  [36]  which  is  typically  much  greater  than  0k  mx.  In  the 

start-up  regime  the  propagation  angle  6k  ^  determines  the  portion  of  the  spontaneous 

power  that  is  within  the  gain  spectrum  and  amplified  as  shown  schematically  in  Fig. 5. 
III.  Transition  from  Incoherent  to  Coherent  Radiation 

The  discrete  nature  of  the  electron  beam  interacting  with  the  pump  laser  field 
leads  to  the  generation  of  spontaneous  (incoherent)  radiation  that  can  be  subsequently 
amplified  [21-26].  During  amplification,  however,  there  is  an  increase  in  the  coherence 
of  the  radiation. 

The  wave  equation  governing  the  x-ray  generation  is 
(V2  -  c~2d2  /3/:)E(r,0  =  4/rc  2g(z)d  J(r, l)/dt  +  4 /rg(z)V p(r,t )  where  the  field  and 
current  density  are  E(r,t),  J(r,r)  =  (E(r,t),J(r,t))e±  +  c.c.,  p{r,t)  is  the  charge  density 
and  g(z)  =  1  for  0  <  z  ^  z  defines  the  interaction  region  and  is  zero  otherwise.  The 
driving  current  density  consists  of  a  coherent  and  an  incoherent  (discrete)  component 

J(x, y,z,t)  =  qY,  v1.(0^-?I(0)^(y->’1(0)^(z-Zi(0) 

(17) 

=  +  qf,  v(°\r0,yo,t)S(r-f'0)(t)), 

is  I 

where  ^"’(f),  vj"’(0  are  the  unperturbed  electron  trajectories  given  by 

v<")  =  (-rao/  yo,  v.o),  JfoA(r,/')  is  the  coherent  current  density  which  is  responsible  for 

the  FEL  interaction  and  the  summation  term  in  Eq.(17)  is  responsible  for  spontaneous 


emission.  Substituting  the  current  density,  Eq.(17),  into  the  wave  equation  and  Fourier 
transforming  the  spatial  and  temporal  variables  we  obtain. 


D(k,co)E(k  ,k  ,k,,co)  = 


%7t2 


X  *  V 

'  Azrq'' 
c  j 


a„  i  


V2  Y0 


ze\p(i  K(k.,co)z/2) 


rs\n(K(k„0))z/2)  ^  N 


K(k.,co)z/2 


X  exP (Ck,./). 


(18) 


i  =  l 


where  K(kt,a>)  =  (co-co0)/\zo  -  (kz+k0) ,  X,.,-  =  “M*.,  “  and  the 

charge  density  has  been  neglected  [23].  In  obtaining  Eq.(l  8)  harmonics  are  neglected 
since  we  assumed  «  1 ,  where  !<5vj  =  an  /(2 yo  k0 )  is  the  magnitude  of  the  electron 

transverse  wiggle  motion  in  the  pump  laser. 

The  function  D(k,ry)  is  given  by 


D(k,a>)  =  - 


Dfel(k,0)) 


(kz  -  {6)/c)(  1  +  (a>-am)/{2fwa>tJff 
(k,  -  fc,(fcx, &?))(£.  -  A2(AX,  &?))(£.  -  k3(kL,oj)) 
(k:+kn-(a)-co(l)/\zo)2 


(19) 


where  the  roots  of  the  dispersion  relation  are  denoted  by  kvk2,  k2  and  A,  denotes  the 
growing  root.  Solving  the  wave  equation  for  the  field  associated  with  the  growing  root, 
i.e.,  integrating  around  the  pole  at  kz  =  A,(Ax,ry),  yields 


E(k.,z,a ))  = 


2(tt)u1  c  r„ 


~—ze\p{iK]z/2) 


^sin(/T,  z/2) A 


K.z/2 


x 


J 


V  “i 

G(kvk1,co )  exp(/A,(A±,ry)z)  exp (/£,.,), 


(20) 


1  =  1 


where  =  K ( A, ,  ry) = ( co-  coo )  /  v  -  ( A,  +  k0 ) = [ ry-  (1  ■ -  f  c2  k\  /  0)2RO )]  1(2 c  f  )  .and 


G  ( A, ,  A± ,  ry)  =  (( A,  +  A0  -  {co- co0 )  /  v  )2 )  /( A,  -  k2  ( A± ,  co)){  A,  -  A,  ( kL ,  ry)) . 
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The  intensity  is  given  by  I(r,z,t)  =  (c/2x)(^E(r,z,t)E* (r,z,t)^  where  (  ) 


denotes  an  average  over  electrons.  If  the  electrons  are  initially  randomly  distributed,  we 


=  Nh  and  obtain 


,2  2 


(E(r,z,a>)E\r,z,a> ))  = 

'  7  4  n  c  Y« 

sin(AT,  z/ 2)' 


f  ki  dki 


zexp 


.K,  z 


,Yi 


K.  z/2 


G(kt ,  k± ,  CD)  Ja ( kLr )  exp(/  /c,  ,  ry)  z) 


(21) 


where  we  have  set  z  -  z  .  The  spectral  power,  defined  by 


dP 

d(0 


=  | 


(£_P_ 

dead  Q.k  J 


OO 

=  — Jr</r(E(r,z,fi>)EV,z,ft>)),  (22) 

o 


is  given  by 


r/2  P 
J  ryr/  Eik 


,  A 


=  8v//2/nc2 
exp(/  Kxz!  2)| 


1  +  a, 


x 


sin(AT.  z/2) 


z/2 


(23) 


IG^,  kL,0))\"  |exp(i  k{ (kL, CO)  z)|\ 


where  kLdk  =  (k  I27t)dilk ,  dilk  is  the  differential  solid  angle  associated  with  the 


wave  vector  and  the  relation  Jrdr  J0(k±  r)  J(l(k'1  r)  =  S(k1-k[)/ k±  was  used. 

o 

a)  Incoherent  Radiation 

The  spectral  brightness  in  the  absence  of  the  FEL  interaction  is  the  spontaneous 
(incoherent)  spectral  brightness  and  is  obtained  from  Eq.(23)  by  setting 
kx  ={0)!c){  1  -  c2  k2 !2(02r)  ,  together  with 
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\exp{iK]Z/2f  =  |G(*,^  ,(o) |  =  |exp(/fc,(fci,<y)z)|2  =1.  The  incoherent  spectral 
brightness  is 


d2P„ 


incoh 


y  d(Odilk  j 


=  8  vfa 


r  _  V 


me 


\K  j 


\2 


1  +  a 


o  J 


sin(ff.z/2) 


K.z/2 


(24) 


The  incoherent  power  radiated  by  the  electron  bunch  per  unit  solid  angle  is 


dP, 


incoh 


dih 


=  Jjry 


i2  A  ~\ 
^  ^inenh 

Kda)dQ.k  , 


= 


vmc 

z 

A 

4 

a. .  * 


V1  +  aw 


(25) 


The  incoherent  power  within  the  solid  angle  AQ.incoh  is 


0  +  «J  4  4 


(26) 


where  A &incuh  <  7t  I  fzo  is  the  solid  angle  associated  with  the  incoherent  radiation  and 
Ph  =vyomci / re  is  the  electron  beam  power. 
b)  Coherent  Radiation 

A  small  portion  of  the  spontaneous  radiation  spectrum  overlaps  the  gain  spectrum 
and  is  amplified  as  depicted  in  Fig.  5.  The  coherently  amplified  portion  of  the  spectrum 
is  determined  by  the  relative  line  widths  of  the  spontaneous  and  gain  spectra,  as  well  as 
the  range  of  amplified  transverse  wave  numbers  kL  given  by  Eq.(16).  The  fractional  line 
widths  associated  with  the  coherent  and  incoherent  (spontaneous)  power  spectrum  are 
respectively 

ScoM{z)l<oK={\l2n)(XJL„)jLjl,  (27a) 


and 


Scinch  (z)  /  a>m  =  0  / 2)  (4  /  z)  ■ 


(27b) 
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The  ratio  of  the  linewidths  is  Sd)a)h  /Seoirunh  =  (]/ x)<Jz/ Lgn  which  implies  that  for 
interaction  distances  less  than  ~  10L^ ,  the  coherent  power  spectrum  is  narrower  than  the 
spontaneous  spectrum.  For  z  >  Lgo ,  we  find  that 


|exp(iA,z/2)|2|  sin(A,  z/2)/(A,  z/2)|  =  (2Z^0/ z )2 ,  and  the  coherent  power  spectral 
brightness  in  Eq.  (23)  is 

/  \2/  \2  /  ,  \2 

exp(r  (cU'±)z),  (28) 


dcodQ.k 


=  (32/9  )vrfmc2 


fz] 

2 

(  \ 
ao 

2 

(L  \ 

R" 

J 

vl+fl«y 

l  ^  J 

where  we  used  G(£,  =1/9,  i.e.,  1/9  of  the  incoherent  power  is  available  for  gain. 

Using  the  power  growth  rate  spectrum  in  Eq.(6),  Eq.(28)  can  be  integrated  over 
frequency  to  give 


dP 


coh 


dnk 


=  6.3  vy^mc2 


„  Y( 


J 


i  \2( 


0+a«  J 


*0 


V  4  J 


A  6)]^ 


ex 


p(r„, -)> 


(29) 


where  A co/ coR0  =  (A0/L  )/ In. 

The  coherent  radiation  beam  is  confined  to  a  narrow  forward  cone  with  solid 
angle  AQ.k  =  x where  04iimil  is  given  by  Eq.(16).  The  incoherent  radiation,  on  the 

other  hand,  is  confined  to  a  cone  angle  6mcoh  ~  1/  yin  >  0kmax. 

The  integration  over  solid  angle  in  Eq.(29)  can  be  approximated  by  evaluating  the 
integrand  at  k±  =  0  and  multiplying  by  the  solid  angle  AQ.k .  The  coherent  power  is 
given  by 


PcAz)  =  25  y> 


I 

f  L  \ 

go 

r, \ 

(L  \ 

go 

0  +  «o2)3 

s  ) 

) 

z  J 

Afl*  Phexp(z/ Lgo). 


(30) 
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The  coherent  power  in  Eq.(30)  will  be  compared  with  GENESIS  simulations  in  the  x-ray 
regime. 

For  interaction  distances  much  greater  than  a  Rayleigh  length,  but  shorter  than  the 
saturation  length,  Lsal  >  z»  ZR  the  propagation  angle  is  equal  to  the  diffraction  angle, 

i.e.,  0k  =  0n  and  the  coherent  power  from  Eq.(30)  becomes 


\2 


Pcoh(z  »  ZR)  =0.Q\ 7] yomc: 


J 


f  2  ,  2  '\ 

h  !rs 

f 


%„/,(z)exp(z/^J,  (31) 


where  Scocoh(z) I coR0  =  (Aco/ coR0)jLJz  =  6.9tjJl^Jz  ,  Tj  =  0.023 (XJ Lgo)  is  the 

conversion  efficiency  and  /  is  the  filling  factor.  The  coherent  power  in  this  limit, 
Eq.(31)  is  similar  in  form  to  that  given  in  [21 ,23,24],  The  ratio  of  the  coherent  power  to 
the  incoherent  power  for  the  same  solid  angle  and  for  z  >  Lgo  is 


p  ' 

 coh 


P  , 

V  incoh  Jmi, 


=  1.6 


'L  ^ 


exp (z/Lg)  =  - 


r L  V/2 

exp(z/  L  ) .  (32) 

^  J 


fa>incoh(z) 

c)  Saturation  length  and  Line  Width 

The  saturation  length  for  the  coherent  radiation  can  be  obtained  by  setting  Pcoh  in 
Eq.(30)  equal  to  the  conversion  efficiency  times  the  electron  beam  power,  tj Pb,  where  tj 
is  given  by  Eq.(13).  The  number  of  power  gain  lengths  at  saturation,  Nmt  =  Lsatl L  ,  is 
given  by 

0  +  fl,2)3  ' 


^/2exp (Nsal)  =  9.2x10 


-4 


(x  1 

2 

f  A0) 

lrJ 

(33) 


The  fractional  line  width  associated  with  the  coherent  radiation,  for  k±  ~  0,  at  saturation 


is 


16 


(34) 


&>«*  =  ri/2Affl  =  N-J/2  K 

®R0  ®R0 

There  is  an  additional  contribution  to  the  line  width  due  to  the  finite  transverse 
wavenumber,  kx ,  spectrum  which  is  given  by  6\ ,  as  indicated  in  Eq.(4). 

IV.  Comparison  of  Theory  with  Simulations 

In  this  section  we  compare  the  analytical  result  for  the  coherent  power,  Eq.(30), 
with  the  simulation  results  from  GENESIS  [33].  GENESIS  simulates  the  conventional 
wiggler  based  FEL  amplifier,  including  start-up.  In  using  GENESIS  to  stimulate  the 
laser  pumped  FEL,  the  wiggler  period  in  GENESIS  is  set  equal  to  twice  the  pump  laser 
wavelength  A*  =  A()/2  =  0.5// m  ,  and  the  wiggler  transverse  gradients  are  removed. 

Besides  GENESIS  there  are  other  FEL  simulation  codes  that  can  be  used  to  simulate  the 
FEL  start-up  physics  [37].  Before  discussing  an  example  of  an  x-ray  laser  pumped  FEL 
it  is  useful  to  consider  the  application  of  the  theory  to  a  conventional  FEL  operating  in 
the  x-ray  regime.  For  this  comparison  we  use  the  Linac  Coherent  Light  Source  (LCLS) 
FEL  atSLAC  [1], 
a)  Wiggler  based  X-Ray  FEL 

The  parameters  of  the  LCLS  FEL  operating  at  15  A  are  given  in  Table  II.  In  the 
GENESIS  simulations  we  use  a  circularly  polarized  wiggler  and  a  cold  electron  beam  to 
make  a  comparison  with  theory.  The  wiggler  strength  parameter  is  therefore  smaller  by  a 

factor  of  ypl  than  the  actual  value  used  in  the  original  LCLS  design  [1],  Figure  6  plots 
the  power  as  a  function  of  propagation  distance  within  the  wiggler.  The  blue  curve  is  the 
result  of  a  Genesis  simulation.  The  dashed  curve  and  solid  black  curves  represent  the 
theoretically  calculated  incoherent  and  coherent  power,  respectively.  The  theoretical 
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incoherent  and  coherent  power  are  functions  of  the  maximum  solid  angle  as  indicated  by 
Eqs.  (26)  and  (30).  In  calculating  the  theoretical  incoherent  power,  we  assume  a 
maximum  solid  angle  of  AQ(Wl  =  7rl{3y.t))2 .  For  the  coherent  power,  we  used  the 

maximum  solid  angle  given  by  Eq.  (16)  and  a  gain  length  consistent  with  the  GENESIS 
simulation.  The  red  curve  represents  the  total  theoretical  power  at  a  given  interaction 
length.  There  is  excellent  agreement  between  theory  and  simulation  in  both  the 
incoherent  (z  <  7  m)  and  coherent  (z  >  7  m)  regimes. 

Figure  7a  plots  the  radiation  intensity  profile  at  saturation  and  shows  that  it  is 
highly  localized  to  region  of  the  electron  beam.  Figure  7b  shows  the  distribution  of  power 
over  transverse  wave-number,  i.e.,  dP / dkx .  Consistent  with  Eq.  (16),  most  of  the  power 
is  contained  with  the  angle  0k>nm ,  which  in  this  case  is  ~  20  times  larger  than  the 

diffraction  angle. 

b)  Laser  Pumped  X-Ray  FEL 

The  parameters  used  in  the  laser  pumped  FEL  are  listed  in  Table  I.  For  these 
parameters  lethargy  effects  (electron  beam  slippage)  are  negligible.  Figure  8  shows  the 
evolution  of  the  x-ray  power  as  a  function  of  interaction  length.  The  theoretically 
calculated  incoherent  and  coherent  powers  are  shown  separately.  In  this  parameter 
regime,  the  maximum  transverse  resolution  of  the  simulation  is  not  sufficient  to  resolve 
the  maximum  propagation  angle  given  by  Eq.  (16).  Hence,  in  calculating  the  theoretical 
incoherent  and  coherent  power  for  comparison  with  the  simulations,  the  maximum  angle 
is  taken  to  be  dsim  -  /l/(2Ax),  i.e.,  the  maximum  angular  resolution  of  the  simulation, 

where  A^  is  the  transverse  grid  size.  The  transition  from  incoherent  to  coherent  radiation 
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occurs  after  ~  2-3  power  gain  lengths.  The  power  gain  length  is  LRo  =  500pm  while  the 
saturation  length  is  Lsal  =8 L  =  0.4cm .  The  conversion  efficiency  is  77  =  0.01  %  which 
corresponds  to  a  saturated  coherent  x-ray  power  of  Pcoh  =  340 kW  .  The  theoretical 
conversion  efficiency,  in  Eq.(13),  gives  a  value  of  0.01%  in  excellent  agreement  with  the 
GENESIS  simulations. 

Figure  9a  plots  the  transverse  x-ray  intensity  profile  at  saturation  and  shows  that 
the  radiation  is  highly  localized  to  the  electron  beam  and  contains  a  large  number  of 
higher  order  transverse  modes.  Figure  9b  plots  the  distribution  of  power  over  transverse 
wavenumber  and  shows  there  is  significant  power  over  the  entire  wave-number  range 
resolved  by  the  simulation. 

In  this  example  the  required  relative  electron  beam  energy  spread  is  <  0.01  %  . 
Higher  electron  beam  energy  spreads  would  substantially  reduce  the  lasing  efficiency  and 
limit  the  growth  of  coherent  x-ray  power.  This  example  indicates  the  stringent 
requirements  placed  on  both  the  electron  beam  and  pump  laser. 

V.  Conclusions 

We  have  analyzed  a  high-gain,  laser  pumped,  x-ray  FEL  amplifier.  The  analysis 
includes  i)  electron  beam  thermal  effects,  ii)  off-axis  propagation  and  iii)  the  transition 
from  incoherent  to  coherent  x-rays.  The  power  gain  length,  saturation  length,  line-width 
and  conversion  efficiency  have  been  calculated  for  the  laser  pumped  FEL.  We  find  there 
is  good  agreement  between  our  theoretical  results  and  GENESIS  simulations.  For 
electron  beams  of  sufficiently  high  quality,  with  energies  of  7  MeV  and  currents  of 
500  A  ,  we  find  that  coherent  x-rays  at  1 5  A  can  be  generated  with  power  gain  lengths  of 
500// m,  saturation  lengths  of  0.4cm  and  conversion  efficiencies  of  0.01%.  To  achieve 
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these  values  the  fractional  electron  beam  energy  spread  must  be  <  0.01  %  .  The  pump 
laser  energy  per  pulse  is  ~  500  J.  While  such  a  coherent  x-ray  source  would  have  a 
number  of  attractive  features,  the  requirements  placed  on  both  the  electron  beam  and 
pump  laser  are  challenging. 
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Electron  Beam  Parameters 

Energy 

Eh  =  7MeV  (yo  =  14.5) 

Current 

Ib  =500  A 

Radius 

rh  =  70// m 

Energy  spread  limit 

t±yly<  j)  =  0.01% 

Pump  Laser  Parameters 

Wavelength 

K  = 

Strength 

a„  =  0-5 

Power 

Po  =  17TW 

Energy/pulse 

*— > 
o 

II 

II 

Pulse  duration 

tl=  2  LstU  !c  =  27  psec 

X-Ray  Parameters 

Wavelength 

II 

Spot  size  (at  saturation) 

rs  =  70  ju  m 

Rayleigh  length  (at  saturation) 

ZR  =  10m 

Power  gain  length 

Lso  =  500// m 

Conversion  efficiency 

77  =  0.01% 

Saturation  length 

Lml  =  0.4  cm 

Saturated  power 

Pm  =  340  kW 

Table  I.  Parameters  for  a  Laser  Pumped-FEL 
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Electron  Beam  Parameters 

Energy 

Eh  =  4.52GeV  (yo  =8.9xl03) 

Current 

Ih=  3.4  kA 

Radius 

rb  =  80// m 

Energy  spread  limit 

A yly<  77  =  0.08% 

Wiggler  (Circular)  Parameters 

Period 

=  3cm 

Strength 

=  2.62 

X-Ray  Parameters 

Wavelength 

lO 

II 

Spot  size  (at  saturation) 

rs  =  140//m 

Rayleigh  length  (at  saturation) 

ZR  =  42m 

Power  gain  length 

L*o  =2m 

Conversion  efficiency 

ri  =  0.08% 

Saturation  length 

Lsa,  =  28m 

Saturated  power 

Psa,  =  12.5  GW 

Table  II.  Parameters  for  LCLS  FEL.  Note  that  the  actual  LCLS  wiggler  is  linearly 
polarized  with  aw  =  3.7  . 
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pump 

laser 


Figure  1:  Schematic  of  laser-pumped  free-electron  laser.  The  pump  laser  and  electron 
beam  propagate  in  opposite  directions  along  the  z  axis. 
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Beam  Current,  Ih  [kA] 

Figure  2:  X-ray  ( A=]5  A)  power  gain  length  Lgo  versus  beam  current  for  a  cold  electron 
beam  with  electron  beam  radius  rfc=50pm  (red)  and  rh  =  75pm  (blue).  The  curves  are 
from  Eq.  (6)  with  Lgo  =  1/F,,,  and  the  solid  circles  are  from  GENESIS  simulations.  The 
other  parameters  are,  a()  =0.5,  =1  pm  and  Eb =7  MeV  . 
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Figure  3:  (a)  Surface  plot  of  normalized  x-ray  growth  rate  T  /r?0  from  Eq.  (8)  versus 
fractional  electron  beam  energy  spread  6ylyo  and  fractional  detuning  (CD -  a)R0) / coR0  for 
the  parameters  of  the  laser-pumped  FEL  of  Table  I.  (b)  Line  plots  of  Ts  /  versus 
Sy/yo  for  various  values  of  detuning  for  the  same  data  as  shown  in  Fig.  3a. 
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Figure  4:  X-ray  conversion  efficiency  versus  beam  current  for  a  cold  electron  beam,  with 
electron  beam  radius  rfr=50pm  (red)  and  rfc=75pm  (blue).  The  curves  are  from 

Eq.(13)  and  the  solid  circles  are  from  GENESIS  simulations.  The  parameters  are  the 
same  as  in  Fig.  2. 
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incoherent  radiation 


Figure  5:  Schematic  diagram  of  the  incoherent  (spontaneous)  and  coherent  (growth  rate) 
spectrum  in  the  ( co ,  k^/k^  plane  showing  the  region  of  overlap.  The  red  area 
corresponds  to  the  region  of  incoherent  emission  of  radiation  from  the  electron  beam 
interacting  with  the  pump  laser.  The  green  area  indicates  the  growth  rate  spectral  region 
in  which  the  radiation  grows  exponentially. 
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Figure  6:  Power  versus  interaction  length  for  the  LCLS  FEL  (cold  beam)  with  parameters 
listed  in  Table  II.  Curves  denote  coherent  power  (solid  black),  incoherent  power 
(dashed),  total  theoretical  power  (red),  and  the  result  of  a  Genesis  simulation  (blue).  The 
theoretical  efficiency  is  T]  =  0.023 (A0  /  Lgo )  =  0.1%.  The  efficiency  observed  in  the 
GENESIS  simulation  is  0.08%. 
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Figure  7:  (a)  GENESIS  simulation  result  showing  transverse  profile  of  intensity  (solid 
curve)  at  z  =  28  m  (saturation)  for  the  LCLS  FEL.  Dashed  curve  denotes  electron  beam 
profile,  (b)  Distribution  of  power  over  normalized  transverse  wavenumber  (propagation 
angle),  for  the  intensity  profile  shown  in  (a).  Propagation  angle  is  normalized  to  the 
diffraction  angle  0D .  0k  tmx  denotes  the  maximum  propagation  gain  angle  from  Eq.  (16). 
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Figure  8:  Power  versus  interaction  length  for  a  laser-pumped  FEL  (cold  beam)  with 
parameters  listed  in  Table  I.  Curves  denote  coherent  power  (solid  black),  incoherent 
power  (dashed),  total  theoretical  power  (red),  and  the  result  of  a  Genesis  simulation 
(blue).  The  theoretical  efficiency  is  r\  =  0.023 (Za/Lgo)  =  0.007%.  The  efficiency 

observed  in  the  GENESIS  simulation  is  0.01%.  In  calculating  the  incoherent  and 
coherent  power  from  Eqs.  (26)  and  (30),  we  used  the  maximum  angle  resolved  by  the 
simulation,  i.e.,  6sim  =  Xl{ 2AX)  =  5x  10  4rad  where  Ax  =  1 .6  pm  is  the  transverse  grid 
size.  In  this  parameter  regime,  0sim  «  6k  mx ,  where  9k  ma  is  the  maximum  propagation 
angle  of  the  coherent  radiation  given  by  Eq.  (16). 
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Figure  9:  (a)  GENESIS  simulation  result  showing  transverse  profile  of  intensity  (solid 
curve)  at  z  =  0.4  cm  (saturation)  for  the  laser-pumped  FEL.  Dashed  curve  denotes 
electron  beam  profile,  (b)  Distribution  of  power  over  normalized  transverse  wavenumber 
(propagation  angle),  for  the  intensity  profile  shown  in  (a).  Propagation  angle  is 
normalized  to  the  diffraction  angle  0D .  Maximum  angle  resolved  by  the  simulation  is 
0k  tm%  =10&D.  Maximum  theoretical  propagation  angle  is  0k  Tmx  ~  500  0P  . 
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